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ABSTRACT 


The systems V949 Cen, V358 Pup, and V1055 Sco are triples comprised of an eclipsing bi¬ 
nary orbiting with a distant visual component on a much longer orbit. The first detailed photometric 
analysis of these interesting systems was performed using also the archival data from Hipparcos, 
ASAS, SuperWASP, OMC, and Pi Of The Sky surveys. The system V358 Pup was also analysed 
using the archival ESO spectra and the radial velocities were derived. The analyses of their light 
curves revealed the physical properties of the eclipsing components, while the interferometric data 
for these systems obtained during the last century show that the binaries are also weakly gravita¬ 
tionally bounded with the third components on much longer orbits. The photometry was carried out 
with the robotic telescope FRAM (part of the Pierre Auger Cosmic Ray Observatory), located in 
Argentina. The BVRI light curves were analysed with the PHOEBE program, yielding the basic 
physical parameters of the systems and their orbits. V949 Cen and V358 Pup were found to be de¬ 
tached systems, while V1055 Sco is probably a semi-detached one. V358 Pup shows a slow apsidal 
motion, while for V1055 Sco we detected some period variation probably due to the third body in the 
system, which cannot easily be attributed to the close visual companion. Therefore, we speculate that 
V1055 Sco can be a quadruple system. For V949 Cen a new orbit was computed, having the orbital 
period of about 855 yr. 

Key words: 

binaries: eclipsing - binaries: visual - stars: fundamental parameters - stars: 
individual: V949 Cen, V358 Pup, V1055 Sco. 

1. Introduction 

The eclipsing binaries as parts of the multiple stellar systems (i.e. of multiplic¬ 
ity three and higher) are excellent objects to be studied. The importance of such 

‘Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under 
programmes ID 082.D-0499, 084.D-0591, 085.C-0614, 086.D-0078, and 190.D-0237. 
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systems lies in the fact that one can study the stellar evolution in them, their ori¬ 
gin, tidal interaction between the components, testing the influence of the distant 
companions to the close pair, a role of Kozai cycles, studying the dynamical effects 
and precession of the orbits, but also the statistics and relative frequency of such 
systems among the stars in our Galaxy (and outside), see e.g. Tokovinin (2007), 
Guinan & Engle (2006), or Goodwin & Kroupa (2005). On the other hand, there 
are still plenty of interesting multiple systems, which were not analysed yet. 

2. Photometric observations 

All of the new photometric observations were obtained in Argentina, using 
the FRAM telescope (Prouza et al. 2010). It is located near the Malargiie town, 
about 1400 m above a sea level. The instrument used for observing V949 Cen and 
V358 Pup was a classical telescope of Schmidt-Cassegrain type, 0.3 m in diameter, 
equipped with a G2-1600 CCD camera. On the other hand, for V1055 Sco only 
a small Ni kk or lens with 107 mm diameter and a CCD camera of G4-16000 type 
was used. Both instruments are equipped with standard Johnson BVRI filters, see 
Bessell (1990). All frames have been reduced by differential aperture photometry 
using software package C-MuniPackQ, the adaptation of the MuniPack code (Hroch 
1998), based on the DaoPhot routines (Stetson 1987). The G2-1600 CCD camera 
sensor area is of 1536x 1024 square pixels (pixel linear dimension is 9 pm), which 
leads to the field of view of 24' xl6' (0.93 arcsec/px), while the second CCD 
camera G4-16000 has the sensor area of 4094x4094 square pixels (pixel linear di¬ 
mension is 9 pm), which leads to the field of view of 7° x7° (6.19 arcsec/px). The 
radii of the apertures for our studied stars were 7.09 px for V1055 Sco, 7.09 px for 
V358 Pup, and 11.82 px for V949 Cen, respectively. The star HD 113916 was used 
as a comparison for V949 Cen, TYC 7636-2147-1 was used for V358 Pup, while 
HD 147704 was used for V1055 Sco. 

The older archival data were taken from the publicly available databases such 
as Hipparcos, ASAS, SuperWASP, OMC, and Pi Of The Sky surveys. However, 
these photometric observations were used mainly for the derivation of the minima 
times for a subsequent period analysis. For the PHOEBE light curve fitting the 
individual filter passbands used were: V for ASAS, H p for Hipparcos, and H p 
also used for SuperWASP due to unavailable SuperWASP filter in PHOEBE and 
relative similarity of transmission curves of H p and SuperWASP filters. 

3. Analysis 

After the reduction of the photometric frames, the FRAM photometry was used 
for a subsequent analysis of the light curves. Program PHOEBE ver. 0.31 (Prsa 
& Zwitter 2005) was used, which is based on a standard Wilson-Devinney code 


*http://c-munipack. sourceforge.net/ 
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Table 1: The individual astrometric observations of V949 Cen. 


Year 

Pos.Angle 

[degree] 

Separation 

[arcsec] 

Reference 

1927.48 

284.9 

0.35 

Finsen (1928) 

1935.84 

275.1 

0.35 

Rossiter (1955) 

1953.37 

277.2 

0.32 

van den Bos (1956) 

1960.49 

267.0 

0.32 

van den Bos (1961) 

1976.128 

248.1 

0.27 

Worley (1978) 

1979.19 

246.0 

0.25 

Heintz (1980) 

1991.25 

231.0 

0.263 

Perryman et al. (1997) 

2011.0401 

206.3 

0.2532 

Hartkopf et al. (2012) 


(Wilson & Devinney 1971) and its later modifications. Because of having only 
limited information about the systems, several parameters of the light curve model 
have to be fixed for the modelling. For the systems with no radial velocities, the 
individual masses of the components are not known. Hence, the mass ratio was 
kept fixed at a value of q = 1.0 for the detached systems (see below). The primary 
temperature was fixed at a value derived from the photometric indices, or according 
to the published spectral type. The albedo and gravity darkening coefficients were 
kept fixed. On the other hand, the value of the third light was also computed due to 
the fact that all of these systems are triple stars. 

4. Individual systems 

For our analysis, we have selected three interesting systems. These systems 
were chosen according to the following selection criteria. All of them are relatively 
bright (<10mag), are on the southern sky, all have periods shorter than 10 days, 
all are multiples (eclipsing binary with a third visual component) and all are also 
only seldom-investigated systems. Therefore, three stars were chosen: V949 Cen, 
V358 Pup, and V1055 Sco. In the following subsections we will focus on the 
individual systems in more detail. 

4.1. V949 Cen 

V949 Cen (also HD 113840, HIP 64025 ) is a variable star of Algol type. Its 
orbital period is of about 3.8 days, while its brightness in maximum is of about 
9 mag in V filter. It was also discovered as a double star by Finsen (1928), and 
since then 8 astrometric observations of the visual pair were carried out, see Table 
[U Houk & Cowley (1975) derived its spectral type as G1/G2V. Its parallax is rather 
uncertain, because the Hipparcos catalogue (van Leeuwen 2007) gave the value of 
parallax 1.42 ±1.48 mas. 

No radial velocity analysis as well as no light curve analysis were published. 
On the other hand, from the available astrometric observations a preliminary visual 
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Figure 1: The orbit of the visual pair of V949 Cen. The individual observations are 
plotted as dots, connected with their theoretical positions on the orbit with short 
abscissae. The dotted line represents the line of the apsides, while the dashed one 
stands for the line of the nodes. The eclipsing binary is located in the coordinates 
( 0 , 0 ). 
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Table 2: The orbital parameters of V949 Cen. 


Parameter 

Value 

Error 

P3 [yr] 

855.3 

36.0 

To [yr] 

1915.4 

36.3 

e 

0.947 

0.006 

a [arcsec] 

2.463 

0.012 

i [deg] 

95.4 

21.3 

Q [deg] 

114.8 

47.6 

oj [deg] 

106.6 

17.5 


orbit was derived, see Hartkopf & Mason (2011). Their analysis resulted in orbital 
period of about 348 yr, but the individual observations rather deviate from the tit 
and the most recent one is very distant from its predicted position. Hence, we 
performed an updated orbital fitting, resulting in parameters given in Table |2] The 
tit is presented in Fig. [Tjand as one can see, it is still very preliminary yet, because 
the arc of the orbit covered with the observations is still rather limited. However, we 
have to mention that the solution is not self-consistent due to high value of the total 
mass from the Kepler’s third law (assuming 3 G1V stars, i.e. 3 • 1.15 M 0 = 3.45 
M 0 ). Hence, the explanation is that the period pj, should be much longer (up to 
the 39000 yr), or other solution is that its distance is much different (i.e. system 
closer, parallax higher). The problematic derivation of parallax from the Hipparcos 
data for the visual binaries was discussed elsewhere (e.g. Soderhjelm 1999), but 
we cannot believe that the difference is so large. 

Concerning the photometry of V949 Cen, we collected all available photome¬ 
try from the Hipparcos (Perryman et al. 1997) and AS AS (Pojmanski 2002) data 
together with our new observations from the FRAM telescopes in BVRI filters. 
There is clearly seen that the system is a detached one of Algol type, having two 
rather similar components. The simultaneous analysis of all available photometry 
was carried out with the PHOEBE code, resulting in light curve parameters given 
in Table [3] and the light curve fits presented in Figure [2] As one can see, the value 
of the third light resulted in a non-negligible contribution, indicating that the third 
component is probably very similar to the eclipsing ones. This is in good agree¬ 
ment with the magnitude difference as presented in already published papers (e.g. 
the Hipparcos catalogue presented the magnitude difference about 0.82 mag). 

The AS AS, Hipparcos and FRAM photometry were also used to derive the 
times of minima for V949 Cen. For computing the times of minima a classical 
Kwee-van Woerden (1956) method was used together with our new AFP method 
(see Zasche et al. 2014). However, the twenty minima times plotted in the O -C 
diagram (see Fig. [3]) show no period variation and the ephemerides suitable for any 
future observations are given in Table [3j Also from this plot it is obvious that the 
innner orbit of the eclipsing pair is perfectly circular. All of the data points used 
for the analysis are stored in the online-only Appendix tables. 
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Phase 

Figure 2: The available light curves of V949 Cen. The ASAS and Hipparcos curves 
were shifted in the y-axis for a better clarity. The bottom plots show the residuals 
of the fit. 
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Figure 3: O-C diagram of V949 Cen. Black symbols stand for the primary min¬ 
ima, the blue ones for the secondary minima. The larger the symbol, the higher the 
weight (higher the precision). 


Table 3: The light-curve parameters of V949 Cen, as derived from our analysis. 


Parameter 

Value 

Parameter 

Value 

HJDq 

2448504.0866 ± 0.0096 

Li, b [%] 

39.9 ± 1.4 

P 1 d | 

3.7867969 ± 0.0000061 

L 2 .b [%] 

30.9 ± 1.1 

i [deg] 

83.66 ± 0.32 

L 3 ,b [%] 

29.2 ± 0.9 

T i [K] 

5900 (fixed) 

U ,v [%] 

38.8 ± 0.9 

T 2 [K] 

5736 ± 42 

L 2 , v [%] 

31.5 ±0.7 


7.429 ±0.221 

L xv [%] 

29.7 ± 0.5 

0.2 

7.679 ± 0.204 

L U r [%] 

38.1 ±0.9 

Fi = F 2 

0 (fixed) 

L 2 ,r [%] 

31.6 ±0.8 

Ai =A 2 

0.5 (fixed) 

L 3 ,r [%] 

30.3 ± 0.5 

Gi = G 2 

0.32 (fixed) 

L u [%] 

37.3 ± 0.7 

q = M 2 /M\ 

1.0 (fixed) 

L 2 j [%] 

31.5 ±0.5 

R\/a 

0.156 ±0.005 

C 3 ,7 [%] 

31.2 ±0.5 

R 2 /a 

0.150 ±0.005 

Ll,ASAS [%] 

38.7 ± 0.7 



R 2 ,asas [%] 

31.6 ±0.6 



L 3 ,asas l%] 

29.7 ± 0.6 



L\jhp [%] 

39.1 ± 1.7 



L 2 jup [%] 

31.4 ±0.9 



L 3 J!IP [%] 

29.5 ± 2.1 


O-C (Period) 
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Phase 

Figure 4: Light curves of V358 Pup. The AS AS and Hipparcos curves were shifted 
in the y-axis for a better clarity. The bottom plots show the residuals of the fit. 

4.2. V358 Pup 

V358 Pup (also FID 51994, F1IP 33487) is an Algol-type eclipsing binary dis¬ 
covered by the Flipparcos satellite (Perryman et al. 1997). Its combined spectral 
type was derived as G5V (Houk 1978). It is also known as a visual double and only 
eight observations were obtained since its discovery in 1926. However, no visible 
motion on the sky is presented, hence its period is probably rather long. The latest 
observation of the double from 1993 indicates the separation of about 1" at a po¬ 
sition angle of 149°. The only detailed study of the star is that one by Clausen et 
al. (2001). However, their uvby light curves were not analysed, only a brief state¬ 
ment about a large fraction of the third light and several times of minima (showing 
mildly eccentric orbit) were published. 

Our new FRAM light curves were analysed simultaneously with the other avail¬ 
able photometry using the PHOEBE code. The results are plotted in Fig. @1 while 
the parameters are given in Table [4| As one can see, the amount of the third light 
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Table 4: The light-curve parameters of V358 Pup, as derived from our analysis. 


Parameter 

Value 

Parameter 

Value 

T i [K] 

5780 (hxed) 

i [deg] 

89.35 ±0.17 

t 2 [ K] 

5157 ±53 

HJD 0 

2450818.7131 ±0.0009 

£41 

22.11 ±0.50 

P 1 d | 

6.7939106 ±0.0000019 

q 2 

20.17 ±0.43 

e 

0.008 ± 0.005 

U ,B [%] 

33.8 ± 1.6 

a>o [deg] 

97.4 ± 3.9 

L 2 ,b [%] 

15.2 ± 1.4 

dm/d t [deg/cycle] 

0.002 ± 0.001 

L xb [%] 

51.0 ±0.5 

F\ = Pi 

0(hxed) 

Li, v [%] 

32.9 ± 1.0 

Ai =A 2 

0.5 (hxed) 

L 2 ,v [%] 

17.7 ± 1.1 

Gi=G 2 

0.32 (hxed) 

L xv [%] 

49.4 ± 0.5 

A |R 0 | 

18.67 ±0.05 

Li,r [%] 

34.3 ± 1.4 

q = M 2 /M\ 

0.88 ± 0.05 

L 2 ,r [%] 

20.4 ± 0.9 

v y [km/s] 

-2.23 ± 0.03 

L X r [%] 

45.2 ± 0.6 

Derived physical quantities: 

L u [%] 

34.5 ± 0.9 

Mi [M q ] 

1.01 ±0.04 

L 2 j [%] 

21.3 ± 1.1 

Ml [Mq] 

0.89 ± 0.03 

L 3 ,7 [%] 

45.2 ± 0.8 

Ri [R q ] 

0.88 ± 0.03 

^l,ASAS [%] 

33.6 ± 2.7 

R 2 [R q ] 

0.86 ± 0.02 

L 2 ,asas [%] 

20.1 ±2.6 

K\ [km/s] 

65.2 ± 2.6 

^3,ASAS [%] 

46.3 ± 4.0 

Ki [km/s] 

73.9 ± 2.7 

L\jup [%] 

32.1 ± 1.1 

M b oi i [mag] 

5.02 ± 0.08 

F 2 ,hip [%] 

17.9 ± 1.2 

M b oi 2 [mag] 

5.56 ± 0.06 

[%] 

50.0 ± 3.3 




Table 5: The individual radial velocities of V358 Pup. 


HJD 

RV i 
[km/s] 

rv 2 

[km/s] 

2454889.75029 

-65.33 

±0.47 

69.54 

± 

2.48 

2454890.72131 

-51.67 

±0.30 

53.90 

± 

0.91 

2455118.89667 

18.27 

±0.44 

-26.06 

± 

0.77 

2455145.81532 

32.13 

±0.19 

-40.86 

± 

0.37 

2455469.82594 

35.89 

± 1.17 

-46.11 

± 

2.81 

2455469.88360 

38.19 

±0.49 

-46.99 

± 

1.51 

2455479.74664 

-26.06 

±0.27 

24.21 

± 

0.48 

2455504.87144 

61.51 

±0.21 

-74.38 

± 

0.34 

2456214.88501 

-66.73 

±0.47 

70.46 

± 

1.42 


is really large, which is in agreement with the published values of the magnitude 
difference between the two visual components (the third component is only about 
0.1 mag fainter than the combined magnitude of the eclipsing pair). 
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Figure 5: The radial velocities of V358 Pup, the black dots and the solid line stand 
for the primary, while the red dots and the dashed line stand for the secondary. 

For the whole analysis we also used the ESO spectra downloaded from the 
online ESO archive (altogether there were used nine HARPS spectra for the anal- 
ysis0). The original data were reduced using the standard ESO routines. The final 
radial velocities (hereafter RV) used for the analysis were derived via a manual 
cross-correlation technique (i.e. the direct and flipped profile of the spectral lines 
manually shifted on the computer screen to achieve the best match) using program 
SPEFO (Horn et al. 1996, Skoda 1996) on several absorbtion lines in the measured 
spectral region (usually Fe and Ca lines). The derived radial velocities are given 
in Table [5] There were also visible some lines from the third body in the spec¬ 
tra, but all these lines show similar radial velocities, of about -2.7 km-s -1 with no 
indication of change during the observed period (of about 3.6 yr). 

As one can see from our plot given in Fig. [51 the final RV curve is still not per¬ 
fect (as it is obvious from the residuals). We suppose that this effect can be partly 
caused by the program PHOEBE itself (and the relative weightening between the 
photometry and spectroscopy) and partly also due to its rather uncertain eccen¬ 
tricity. From the pure RVs it would indicate that the orbit has lower eccentricity, 
but from the photometry and the minima times there resulted that the eccentricity 
should be higher. 

As one can see from the final parameters given in Table [4] the errors of masses 

*HARPS is a spectrograph mounted on the 3.6-m ESO telescope at the La Silla Observatory. Its 
spectral range is from 3781 to 6911 A and the resolving power is 115 000. 
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Figure 6: Evolutionary tracks for the system V358 Pup and both its components, 
star 1 (primary in black) and star 2 (secondary in red) plotted for three different 
metallicities. The blue dashed line represents the ZAMS line for Z = 0.02. 


2448000 2450000 2452000 2454000 2456000 2458000 2460000 



Figure 7: O — C diagram of V358 Pup, showing very slow apsidal motion. Black 
symbols stand for the primary minima, the blue ones for the secondary minima, 
and the larger symbols correspond to the data points with the higher weights. 
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and radii are very low and one has to doubt about their values. With only a few data 
points and our poor fit one cannot probably derive a mass with a precision of about 
3-4%. Therefore, we would like to point out that the presented solution provided 
with the PHOEBE program is a formal one, based on the abovementioned assump¬ 
tions (the inclusion of some future knowledge of these stars can significantly shift 
our solution). Moreover, also the errors as resulted from the code are purely math¬ 
ematical ones and usually are strongly underestimated in the PHOEBE program 
(Prsa & Zwitter 2005, or PHOEBE manuad). 

Another indication that the presented errors of parameters as given in Table 0] 
are underestimated is a plot in Fig. [6] The evolutionary tracks are plotted together 
with the position of the two eclipsing components as resulted from our analysis 
(values of temperatures and bolometric magnitudes taken from Table 0] with the 
error of T\ equal to error of 73) and the theoretical curves are taken from the EZ- 
wetH which is using a well-known MESA code (Paxton et al. 2011). To conclude, 
here we deal very probably with two stars still located on the main sequence or its 
vicinity and the models indicate that the stars have a near solar metallicity. 

From the available photometry also the minima times were derived. We col¬ 
lected the most recent FRAM photometry together with the already published min¬ 
ima from the study by Clausen et al. (2001). Moreover, also the older Hipparcos 
(Perryman et al. 1997) and ASAS (Pojmanski 2002) data were used. All of these 
photometric data yielded 26 minima times, however some of them (especially the 
ones from the Hipparcos and ASAS data) are of poor quality only. The resulting 
O-C diagram is plotted in Figure [7] The ephemerides given in Table 0]refer to the 
y = 0 line in Fig. 0) We analysed these data using the apsidal motion hypothesis 
as described e.g. by Gimenez & Garcfa-Pelayo (1983). This led to the finding that 
the apsidal motion is very slow, having the period of thousands of years and our 
present analysis is still only very preliminary yet (and as one can see, the errors 
of apsidal motion parameters are still very large). From our apsidal motion analy¬ 
sis there resulted that the relativistic contribution to the total apsidal motion rate is 
rather small, being only about 11 %, and the value of log/G = -0.165. This value 
can be compared with the evolutionary tracks as published e.g. by Claret (2004), 
resulting in a finding that such a log A :2 value can be achieved with a 1 M 0 star with 
an age of about 2.5 • 10 9 yr. 

4.3. VI055 Sco 

The system V1055 Sco (also HD 148121, HIP 80603) is the only one W UMa- 
type binary in our sample of stars. It has the orbital period of about 0.36 days only 
and rather shallow amplitude of photometric variation of about 0.25 mag. It was 
discovered from the Hipparcos data (Perryman et al. 1997). On the other hand, 
Szczygiel et al. (2008) included the star into their compilation of the ASAS eclips- 


■''http://phoebe-project.org/1.0/docs/phoebe _ manual.pdf 
® http://www.astro.wisc.edu/ ~ townsend/static. php?ref=ez-web 
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Table 6: The light-curve parameters of V1055 Sco, as derived from our analysis. 


Photometry 

FRAM 

AS AS 

SuperWASP 

Hip 

Parameter 

Value 

Value 

Value 

Value 

HJDq 


2452050.7441 ±0.0043 


P 1 d | 


0.3636835 

± 0.0000038 


T i [K] 


5850 (fixed) 


T 2 [K] 

5751 ±30 

5632 ± 26 

5881 ±23 

5499 ± 75 

i [deg] 

89.4 ±0.8 

88.2 ± 1.3 

88.6 ±0.8 

87.0 ± 1.9 


6.391 ±0.019 

6.746 ±0.032 5.917 ±0.021 

5.407 ±0.105 

Q.2 

6.293 

6.422 

5.782 

4.410 

q = M 2 /M\ 

3.0 ±0.1 

3.1 ±0.1 

2.7 ±0.1 

1.7 ±0.1 

U [%] 

13.5 ± 0.2 (B) 

12.8 ±0.5 

12.5 ± 0.3 

14.3 ±0.4 

U [%] 

13.2 ± 0.2 (V) 




U [%] 

13.0 ± 0.2 (R) 




U [%] 

12.9 ± 0.2 (I) 




l 2 [%] 

37.1 ±0.2 (B) 

37.8 ± 0.5 

39.3 ± 0.4 

35.0 ±0.8 

l 2 [%] 

37.3 ± 0.2 (V) 




l 2 [%] 

37.4 ± 0.2 (R) 




l 2 [%] 

37.4 ± 0.2 (I) 




Li [%] 

49.4 ± 0.3 (B) 

49.4 ± 0.7 

48.3 ± 0.6 

50.7 ± 0.9 

L 3 [%] 

49.5 ± 0.3 (V) 




Li [%] 

49.6 ± 0.2 (R) 




Li [%] 

49.7 ± 0.2 (I) 




F i =F 2 


0.0 

(fixed) 


Ai=A 2 


0.0 

(fixed) 


G\ = G 2 


1.0 

(fixed) 


R\ja 

0.297 ± 0.022 

0.277 ± 0.017 0.303 ± 0.012 

0.263 ± 0.019 

R 2 /a 

0.517 ±0.045 

0.517 ±0.035 0.502 ±0.029 

0.452 ± 0.040 


ing binaries with the coronal activity. Torres et al. (2006) presented the spectral 
type of V1055 Sco as G2V. The first astrometric observation of the visual double 
was carried out in 1925, published by Finsen (1941). Since then 13 positional mea¬ 
surements were obtained, but the movement on the orbit is very slow (orbital period 
of the order of thousands of years probably). The latest observation of the pair from 
2010 indicates that the separation is of about 0.125" at a position angle of 217°. 
Magnitude difference between the two visual components was derived to be from 
0.0 to 0.3 mag (slightly different in different pass bands). 

The light curve solution was a bit more complicated here than for the previous 
two systems. At first, even the ephemerides were not very certain due to the fact 
that in our new FRAM photometry both eclipses seem to be of similar depths (i.e. 
both temperatures probably very similar), but the Hipparcos data show much larger 
difference between the minima depths. Therefore, we decided to use the deeper 
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Table 7: The LITE orbital parameters of V1055 Sco. 


Parameter 

Value 

Error 

P3 [yr] 

36.7 

12.5 

A [day] 

0.031 

0.009 

e 

0.00 (fixed) 

f(m 3 ) [M q ] 

0.118 

0.037 


Hipparcos minimum as the primary one and set the linear ephemerides according 
to this. For the analysis and subsequent derivation of minima times for the period 
analysis we collected the Hipparcos, ASAS, SuperWASP (Pollacco et al. 2006), Pi 
of the sky (Burd et al. 2005), OMC (Mas-Hesse et al. 2004), and our new FRAM 
data. From these different data sets the ASAS, SuperWASP, and Hipparcos seem to 
be the best ones for an analysis, however the total eclipses are only hardly visible. 

We solved the individual light curves of V1055 Sco separately. This is due to 
the fact that the system undergoes some changes of the orbital period between the 
different datasets (see below) and for the individual data the different ephemerides 
were used (fixed from our solution from the O-C diagram below). The other two 
systems (V949 Cen and V358 Pup) were analysed with the complete photometry 
simultaneously, because it is much more robust. For the result see Table |6] with 
the final parameters of our fit. As one can see from the table, the results are rather 
inconsistent with each other. The detached or contact configuration of the binary 
did not led to an acceptable result. We applied a classical "q-search method" for 
deriving the mass ratio of the binary (this can be applied because we deal with a 
semidetached system with secondary filling its respective Roche lobe). The param¬ 
eters as derived from some of the data sets are too different from the others (e.g. 
the temperatures), that one has to doubt about conclusiveness of the results. Gen¬ 
erally, one can conclude that the role of primary and secondary is very probably 
interchanged (due to mass ratio and luminosity), but the more massive component 
is probably the cooler one. Also rather remarkable is quite high contribution of the 
third light in the system. We can speculate about the origin of such different results 
based on different data sets. One possible explanation is that the system is active 
and the shape of the light curve varies in time. A role of some spots on the surface 
is also not excluded (slight asymmetry of the curve). What is also remarkable is the 
fact that the values of inclination seem to be increasing in time, but only a dedicated 
future observation can prove or rule out this hypothesis. 

From the available photometry, we also derived the minima times for a period 
analysis. Plotting the O-C diagram we found that there is a weak sinusoidal 
variation. Therefore, we applied a hypothesis of the third body modulating the 
orbital period of the eclipsing pair, a so-called Light-Time Effect (hereafter LITE), 
see e.g. Mayer (1990) or Irwin (1959). The result is given in Table |7] and the 
final fit is presented in Fig. [9] The same as for V358 Pup applies here and the 
ephemerides given in Table [6] represent the linear ephemerides in Fig. [9] As one 





14 


A. A. 



-0.4 -0.2 0 0.2 0.4 0.6 0.8 

Phase 


Figure 8: Light curves of V1055 Sco. The ASAS and Hipparcos curves were 
shifted in the y-axis for a better clarity. The bottom plots show the residuals of the 
fit. 
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Figure 9: The O-C diagram of V1055 Sco, plotted with all available times of 
minima observations. The dots stand for the primary, while the open circles for the 
secondary minima. The larger symbols stand for the observations with the higher 
precision. The solid curve represents the final LITE fit. The bottom plot shows the 
residuals of the fit. 

can see, the coverage of the orbit is still rather poor, hence we used only a circular 
orbit as a fit. If we collect more observations in the upcoming years, some more 
detailed analysis would be possible. The errors of the parameters are so large due 
to insufficient number of data points, especially close to the maximum in the O-C 
diagram. Hence, the amplitude as well as the orbital period can only be estimated. 
We also tried to identify this period variation with the visual component observed 
via interferometry, but this attempt was not successful at all. Much different periods 
and eccentricities are needed for both orbits to adequately fit the data. This would 
indicate that here we probably deal with a quadruple stellar system. Some future 
observations in upcoming years would help us to identify the period variation with 
the visual component or to prove its quadruple nature. 

5. Discussion and Conclusions 

We performed the first analysis of the light curves for three southern eclips¬ 
ing binaries containing the third components. However, these distant components 
have rather long orbital periods, the only exception is V949 Cen, where we pre¬ 
sented a preliminary orbital solution of the visual double. All of these systems 
were observed remotely with the FRAM telescopes located in Argentina and the 
light curve analysis revealed some interesting results. From the light curve fits we 
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also tried to derive the individual colors for the components and from the photo¬ 
metric indices B-V derive the temperatures. Such values can be compared with 
the temperatures as resulted from the light curve fits to confirm whether the results 
are self-consistent. We found that the method is not sensitive enough and its results 
are influenced by relatively large errors. However, the resulted temperatures agree 
with the ones from the light curve fits given in Tables [3jSJ and[6]within 100-150 K. 

Probably the most problematic issue of our presented study is the fact that for 
the whole analysis we assumed the value of primary temperatures for the light 
curve solution based only on the photometric indices or the published spectral 
types. However, this can be quite problematic due to the fact that for V358 Pup 
and V1055 Sco the most luminous components seem to be their tertiaries. 

V949 Cen and V358 Pup were found to be well-detached systems, while the 
system V1055 Sco is probably a semi-detached one. V358 Pup is showing slow 
apsidal motion of the inner orbit, but its period is very long. For V1055 Sco we 
presented a preliminary solution for the period variation, which can possibly be 
confirmed via a dedicated observation via spectroscopy and/or interferometry. 
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Table 8: Heliocentric minima times for the studied systems. 


Star 

HJD 

2400000+ 

Error 

[days] 

Source 

Star 

HJD 

2400000+ 

Error 

[days] 

Source 

V949Cen 

48492.72269 

0.00564 

HIP 

V1055 

Sco 

53909.70012 

0.00085 

Pi of the sky 

V949 Cen 

48494.62638 

0.00530 

HIP 

V1055 

Sco 

53910.60608 

0.00081 

Pi of the sky 

V949Cen 

52199.99788 

0.00455 

ASAS 

V1055 

Sco 

53938.61031 

0.00141 

Pi of the sky 

V949 Cen 

52201.89344 

0.00458 

ASAS 

V1055 

Sco 

53941.70325 

0.00112 

Pi of the sky 

V949Cen 

52908.13137 

0.00287 

ASAS 

V1055 

Sco 

54257.02189 

0.00084 

ASAS 

V949 Cen 

52910.02381 

0.00170 

ASAS 

V1055 

Sco 

54257.20413 

0.00120 

ASAS 

V949 Cen 

53669.27729 

0.00242 

ASAS 

V1055 

Sco 

54626.89174 

0.00121 

ASAS 

V949Cen 

53671.17083 

0.00712 

ASAS 

V1055 

Sco 

54627.07413 

0.00150 

ASAS 

V949 Cen 

54627.33616 

0.00742 

ASAS 

V1055 

Sco 

54986.40150 

0.00080 

ASAS 

V949Cen 

54629.22726 

0.00306 

ASAS 

V1055 

Sco 

54986.58393 

0.00108 

ASAS 

V949 Cen 

56660.84562 

0.00529 

FRAM 

V1055 

Sco 

56554.62091 

0.00481 

FRAM 

V949Cen 

56677.89022 

0.00629 

FRAM 

V1055 

Sco 

56557.53361 

0.00261 

FRAM 

V949 Cen 

56679.78165 

0.00474 

FRAM 

V1055 

Sco 

56559.53170 

0.00300 

FRAM 

V949Cen 

56681.67097 

0.00867 

FRAM 

V1055 

Sco 

56762.83229 

0.00217 

FRAM 

V949 Cen 

56696.82051 

0.00050 

FRAM 

V1055 

Sco 

56763.74187 

0.00219 

FRAM 

V949 Cen 

56736.58168 

0.00297 

FRAM 

V1055 

Sco 

56781.74747 

0.00271 

FRAM 

V949Cen 

56791.48915 

0.01289 

FRAM 

V1055 

Sco 

56835.75144 

0.00502 

FRAM 

V949 Cen 

56823.68156 

0.01061 

FRAM 

V1055 

Sco 

56836.66608 

0.00317 

FRAM 

V949Cen 

56840.72217 

0.00448 

FRAM 

V1055 

Sco 

53862.41972 

0.00219 

SuperWASP 

V949Cen 

56844.50644 

0.01036 

FRAM 

V1055 

Sco 

53862.59620 

0.00082 

SuperWASP 

V358 Pup 

48440.84746 

0.00202 

HIP 

V1055 

Sco 

53881.33135 

0.00101 

SuperWASP 

V358 Pup 

48444.23326 

0.00163 

HIP 

V1055 

Sco 

53881.50766 

0.00072 

SuperWASP 

V358 Pup 

52517.19479 

0.00046 

ASAS 

V1055 

Sco 

53891.51779 

0.00656 

SuperWASP 

V358 Pup 

52520.58192 

0.00071 

ASAS 

V1055 

Sco 

53909.33457 

0.00315 

SuperWASP 

V358 Pup 

53624.60224 

0.00245 

ASAS 

V1055 

Sco 

53921.33692 

0.06034 

SuperWASP 

V358 Pup 

53627.99432 

0.00671 

ASAS 

V1055 

Sco 

53961.34467 

0.00046 

SuperWASP 

V358 Pup 

54630.10122 

0.00204 

ASAS 

V1055 

Sco 

54207.56117 

0.00161 

SuperWASP 

V358 Pup 

54633.48790 

0.00248 

ASAS 

V1055 

Sco 

54208.46336 

0.00023 

SuperWASP 

V358 Pup 

50784.74575 

0.00030 

Clausen et al. (2001) 

V1055 

Sco 

54211.56150 

0.00068 

SuperWASP 

V358 Pup 

50801.72630 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54212.46443 

0.00041 

SuperWASP 

V358 Pup 

50818.71420 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54212.65260 

0.00105 

SuperWASP 

V358 Pup 

50835.69555 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54216.46515 

0.00054 

SuperWASP 

V358 Pup 

50869.66590 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54216.64935 

0.00240 

SuperWASP 

V358 Pup 

51110.85360 

0.00100 

Clausen et al. (2001) 

V1055 

Sco 

54230.47724 

0.00387 

SuperWASP 

V358 Pup 

51127.83470 

0.00030 

Clausen et al. (2001) 

V1055 

Sco 

54231.56459 

0.00383 

SuperWASP 

V358 Pup 

51134.62900 

0.00100 

Clausen et al. (2001) 

V1055 

Sco 

54236.65049 

0.00081 

SuperWASP 

V358 Pup 

51144.82295 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54237.38584 

0.00056 

SuperWASP 

V358 Pup 

51151.61695 

0.00020 

Clausen et al. (2001) 

V1055 

Sco 

54237.55969 

0.00436 

SuperWASP 

V358 Pup 

51212.76235 

0.00030 

Clausen et al. (2001) 

V1055 

Sco 

54238.47345 

0.00231 

SuperWASP 

V358 Pup 

51236.53670 

0.00100 

Clausen et al. (2001) 

V1055 

Sco 

54239.38623 

0.00183 

SuperWASP 

V358 Pup 

56379.52621 

0.00288 

FRAM 

V1055 

Sco 

54244.29166 

0.00010 

SuperWASP 

V358 Pup 

56661.48240 

0.00482 

FRAM 

V1055 

Sco 

54246.29080 

0.00137 

SuperWASP 

V358 Pup 

56664.86830 

0.00153 

FRAM 

V1055 

Sco 

54246.47502 

0.00750 

SuperWASP 

V358 Pup 

56688.65368 

0.01859 

FRAM 

V1055 

Sco 

54247.37822 

0.00218 

SuperWASP 

V358 Pup 

56722.62160 

0.00345 

FRAM 

V1055 

Sco 

54247.56558 

0.00169 

SuperWASP 

V358 Pup 

56756.59497 

0.05620 

FRAM 

V1055 

Sco 

54248.29360 

0.00082 

SuperWASP 

VI055 Sco 

48048.77191 

0.00165 

HIP 

V1055 

Sco 

54266.30300 

0.00151 

SuperWASP 

VI055 Sco 

48048.95356 

0.00095 

HIP 

V1055 

Sco 

54268.29609 

0.00450 

SuperWASP 

VI055 Sco 

48674.66215 

0.00694 

HIP 

V1055 

Sco 

54268.47414 

0.00071 

SuperWASP 

VI055 Sco 

48674.84790 

0.00250 

HIP 

V1055 

Sco 

54271.20528 

0.00030 

SuperWASP 

VI055 Sco 

52050.71267 

0.00081 

ASAS 

V1055 

Sco 

54271.38183 

0.00242 

SuperWASP 

VI055 Sco 

52050.89202 

0.00089 

ASAS 

V1055 

Sco 

54271.56979 

0.00004 

SuperWASP 

VI055 Sco 

52501.32187 

0.00199 

ASAS 

V1055 

Sco 

54272.29748 

0.00202 

SuperWASP 

VI055 Sco 

52501.49955 

0.00157 

ASAS 

V1055 

Sco 

54272.47287 

0.00102 

SuperWASP 

VI055 Sco 

52795.90580 

0.00071 

ASAS 

V1055 

Sco 

54273.38677 

0.00393 

SuperWASP 

VI055 Sco 

52796.08872 

0.00063 

ASAS 

V1055 

Sco 

54274.29147 

0.00205 

SuperWASP 

VI055 Sco 

53090.1272 

0.0025 

Zasche et al. (2009) 

V1055 

Sco 

54274.47862 

0.00243 

SuperWASP 

VI055 Sco 

53090.3076 

0.0031 

Zasche et al. (2009) 

V1055 

Sco 

54288.47677 

0.00439 

SuperWASP 

VI055 Sco 

53113.04009 

0.00149 

ASAS 

V1055 

Sco 

54294.29489 

0.00034 

SuperWASP 

VI055 Sco 

53113.21796 

0.00163 

ASAS 

V1055 

Sco 

54537.60107 

0.00159 

SuperWASP 

VI055 Sco 

53429.09480 

0.00140 

OMC 

V1055 

Sco 

54562.51977 

0.00294 

SuperWASP 

VI055 Sco 

53483.6384 

0.0013 

Ogloza et al. (2008) 

V1055 

Sco 

54564.51429 

0.00065 

SuperWASP 

VI055 Sco 

53483.8200 

0.0005 

Ogloza et al. (2008) 

V1055 

Sco 

54586.52250 

0.00185 

SuperWASP 

VI055 Sco 

53530.37204 

0.00091 

ASAS 

V1055 

Sco 

54588.34293 

0.00106 

SuperWASP 

VI055 Sco 

53530.55329 

0.00079 

ASAS 

V1055 

Sco 

54594.33636 

0.00549 

SuperWASP 

VI055 Sco 

53869.88015 

0.00288 

ASAS 

V1055 

Sco 

54601.42667 

0.00723 

SuperWASP 

VI055 Sco 

53870.05697 

0.00184 

ASAS 

V1055 

Sco 

54614.52627 

0.00086 

SuperWASP 

VI055 Sco 

53909.51280 

0.00155 

Pi of the sky 











